arXiv:1504.03594vl [stat.ME] 14 Apr 2015 


Splitting hybrid Make-To-Order and Make-To-Stock demand profiles 

Wolfgang Garn,James Aitken 

The Surrey Business School, University of Surrey, Guildford, Surrey, GU2 7XH, United Kingdom 


Abstract 

In this paper a demand time series is analysed to support Make-To-Stoek (MTS) and Make-To-Order 
(MTO) produetion decisions. Using a purely MTS production strategy based on the given demand can lead 
to unnecessarily high inventory levels thus it is necessary to identify likely MTO episodes. 

This research proposes a novel outlier detection algorithm based on special density measures. We divide 
the time series’ histogram into three clusters. One with frequent-low volume covers MTS items whilst a 
second accounts for high volumes which is dedicated to MTO items. The third cluster resides between the 
previous two with its elements being assigned to either the MTO or MTS class. The algorithm can be applied 
to a variety of time series such as stationary and non-stationary ones. 

We use empirical data from manufacturing to study the extent of inventory savings. The percentage of 
MTO items is reflected in the inventory savings which were shown to be an average of 18.1%. 

Keywords: demand analysis; time series; outlier detection; production strategy; Make-To-Order(MTO); 
Make-To-Stock(MTS); 


1. Introduction 


Research into production systems has generally 
characterized and modeled them as either make-to- 
order (MTO) or make-to-stock (MTS). We propose 
the following definitions. A make-to-order product 
is a product where the required quantity is manufac¬ 
tured after receiving a sales order. MTO products 
are commonly identified by low average demand an d 


a high coefficient of variation (ISoman et al.L l2007h . 


A make-to-stock product is a product where items 
are manufactured on anticipated demand. A concrete 
sales order may not exist. Usually the required quan¬ 
tity is derived from forecasts. 

Presenting the operating choice of production sys¬ 
tems as a choice between MTO or MTS simplifies 
the discussion but does not reflect current production 
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dynamics. Fewer and fewer firms can be classified as 
purely MTS or MTO in practice ( Christopher . 2010 : 


Soman et ah . 20061: Aitken et ah . 2003). Investigat¬ 


ing and modeling the combined MTS-MTO chal¬ 
lenge has had limited research with onl y a few papers 


exploring some aspects of the problem (ISoman et al. 


20061: Rajagopalan , 2002h . Where capacity exists for 


manufacturers it can be expedient to separate prod¬ 
ucts between MTO and MTS in terms of planning 
and control. The literature highlights the benefits of 
a focusing/isolating production of MTO and MTS in 
terms of changeover s, bottleneck reduc tion, process 
variation and costs dSchmenneii 120041) . However, 
separating products between the two approaches is 
not possible for many firms. Sharing and schedul¬ 
ing finite capacity between MTS and MTO prod- 
ucts is a major c hallenge for many organisations 
(IKerkkanenL [20071). This paper investigates how or¬ 
ganisations manage the MTS-MTO dynamic in the 
context of a food manufacturing business. The food 
sector predominately services supermarkets and is 
under growing pressure to increase the number of 
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stock keeping units (SKUs), reduee lead-times and 
manage the delivery of an inereasing heterogeneous 
serviee demand aeross eustomers. The growing eom- 
plexity of operating in the food sector has increased 
the ineidenees of shared capaeity between MTO and 
MTS. The seetor offers the opportunity of investigat¬ 
ing the eombined MTS-MTO seheduling challenge 
in the context of finite eapacity faeed with high de¬ 
mand variability produets classified as either MTO 
or MTS. 

Within the food seetor eompanies arrange bulk 
sales, promotions or marketing events that impaet on 
the demand pattern of produets. These events are 
not always visible to the production planners lead¬ 
ing to unplanned surges in demand creating unstable 
produetion schedules and diminishing service levels. 
Unplanned surges can lead to produet being purely 
produeed via a MTS process when they aetually re¬ 
quire a MTO produetion strategy. Ignoring this of¬ 
ten results in the depletion of MTS produeed inven¬ 
tory for regular sales and subsequently unplanned 
ehanges in produetion ( Sehmennen. 20041) . We intro- 
duee a new approaeh, whieh identifies MTO quanti¬ 
ties within an otherwise MTS classified product. The 
objective of our paper is to present a method that 
can support the identification of changes in status 
between MTS and MTO. Identification of changes 
can motivate a MTO strategy for produets which are 
treated as MTS, a hybrid approaeh, providing the ad¬ 
vantage of a produetion sehedule whieh reduees in¬ 
ventory and improves eustomer serviee. The ease 
study demonstrates that applying our new method 
led to redueed inventory levels of 18.1% on aver¬ 
age for previously pure MTS classified produets. The 
proposed new method automatically identifies MTO 
items as “outliers” within a time series, allowing 
the absence of order reeords. The novelty of the 
method is the usage of elustering demand frequency 
extracted from time series. 

Our paper is organized as follows; after intro- 
dueing the ease study firm we review the literature 
on MTS-MTO and MTO/MTS in the seetion three 
and then turn in seetion four to the new proposed 
method for the identifieation of change in produet 
status between MTS, MTO and MTO/MTS. Section 
five presents the results from the modeling and im¬ 
plementation of the new approaeh in the ease study 



Raw Material Cleaning & sorting Intermediary Packaging Finished Goods 


Figure 1; Production process. 

firm followed by discussions on the significanee of 
identifying and managing the status ehange between 
MTO and MTS. Seetion six presents the eonelusions 
and suggestions for future work in this area of grow¬ 
ing importanee for operation managers. 

2. Company and Process Overview 

The ease study firm produees over 250 different 
products derived from five raw materials. The raw 
materials are eleaned and sorted leading to the pro¬ 
duetion of eight intermediary semi-finished products 
which are stored in material handling silos awaiting 
release to the appropriate produetion lines for paek- 
ing. Storage of the semi finished produets was lim¬ 
ited to 48 hours due to baeterial and deeay issues. 
The variety of different paekaging options generated 
in exeess of 250 finished produets for sale to the re¬ 
tail and entering market segments. The production 
process for the food processing business is shown in 
figure 1, The manufaeturing process had been faced 
with the challenges of growing sales, inereased sup¬ 
ply ehain eomplexity and falling serviee levels of a 
broadening heterogeneous eustomer base. The result 
of the inereasingly variable demand on the manufae¬ 
turing operation was the eonversion of the produc¬ 
tion and planning control schedule from a plan with 
a three week time horizon to a reaetive customer or¬ 
der list. The reaetive planning and control modus 
operand! had foreed manufacturing into an increas¬ 
ing number of unplanned ehangeovers and diminish¬ 
ing produetivity, effectively reducing the available 
eapacity in the paekaging operation. The firm had 
moved from a position of flexible spare eapaeity to a 
bottleneek over a three year period due to inereased 
eomplexity and reduetion in sehedule stability. 

The ease study eompany refieets the food man¬ 
ufaeturing seetor whieh has eome under inereasing 
pressure from supermarkets and distributors to di¬ 
minish the homogenisation of their offerings, re¬ 
duee lead-times, increase promotions, shorten prod¬ 
uct life-cycles and provide more frequent deliv- 
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cries (Fisher et al. 


19991: iTaylor and Fearnd. l2006l: 


Squire et all l2009f) . The continually changing envi¬ 


ronment that the food manufacturers operate within 
has made the management of processes through 
the application of forecasting methods difficult. 
The inaccuracies of forecasting in the food sec¬ 
tor have been found to have a significant impact 
on the performance and e fficiency of the op eration 
( Tavlor and Fearne . 2006 : Aghazadeh . 2004 ). Man¬ 
ufacturing a large number of products through a 
pure MTS approach and putting them into inven¬ 
tory is not via ble because of unpr edictable demand 
and shelf-life (ISoman et al.L 120071). Due to increas¬ 
ing volatility the effectiveness and returns from in¬ 
creased investment on fore casting could be viewed 
as futile ( Christopher . 2010l) rendering the pure MTS 
approach ineffective. In order to improve the perfor¬ 
mance of the operation the case study firm decided 
to alter its approach to planning and control of MTO 
and MTS categorised products 

3. Literature Review 

3.1. MTS and MTO 

Make-To-Stock (MTS) systems are described in 
terms of low variety a nd high volume that are 


driven through fore casts (iRafiei and Rabbani 12011 


Soman et al.L 120041) . Operationally the issues that 


are in focus for this system are lot size determi- 
natio n, accuracy o f forec asting and inventory con¬ 


trol ( Soman et ah . 20061) . A MTS approach is 


used for manufacturing items which are standard¬ 
ized with high volume, regular d emand and a re de- 


livered to customers from stock (BirouetaL, 2011 


Chang and Lu . 2010 : Kerkkanenl 2007 ). The pre¬ 
ferred option for firms is to move from a MTS to a 
MTO system of production linked to a c hange in phi¬ 
losop h y from a push to pull approach dBirou et ah . 


H K 

20111: iJodlbauerl 120081) . MTO products are out¬ 


lined in terms of customised ite ms with low ayerage 
volumes and irre gular demand ( Birou et^ 2011 


Kerkkanen . 2007b . The uncertainty linked to MTO 
products places planning and control process in a 
pivotal and time-s ensitive critical role within firms 


( Corti et ah . 2006h . Reducing lead-times and im- 


of demand uncertainty, has been a focus for sev 
eral researchers in developing models that can sup 
port firms in their decision making t o optimize pro 
ducticin schedules and performance (IZaerpour et al. 
2009; .lodlbauer . 2008 : lsa\^ . 2006b . 


An alternative approach to managing the pres¬ 
sure related to reducing lead-times and growing un¬ 
certainty in demand is through the development of 
a hybrid MTS/MTO approach. The approach at¬ 
tempts to harness the strengths of the pure MTO 
and MTS systems (fFedergruen and Katalan . 1999b . 
Through operating a MTS system in the downstream 
sections of the operations parts can be produced 
that are subsequently assemble d in the MTO sec¬ 
tion w hen orders are received (IRafiei and RabbaniL 
201 ill) . Heuristic models have been developed to 


control and schedule production activities to support 
the portioning pro blem and successfu l ly operate the 


hybri d approach (IChang et al.l. l2003l: iRajagopalanl. 


2002b . The splitting of the operation reflects the 
concept of “postponement which has successfully 
been applied in many sectors where it is possible to 
store semi-finished products or modular parts await¬ 
ing orders before assembling the fin al product with¬ 


out any detritus effect on quality (IZaerpour et al. 


20091: iKerkkanenl. l2007h . The option to build semi¬ 


proving the reliability of due dates, in the context 


finished product and store the materials until an order 
is received has limited application in the food pro¬ 
cessing sector. Strict food storage guidelines, bac¬ 
terial control challenges and short shelf-life render 
the opportunity for utilizing a hybrid MTS/MTO ap¬ 
proach, th rough semi-finished p roducts, difficult if 
not illegal ( Johnston et al.l 2003 ). However the mix¬ 
ing of the MTO and MTS approaches has merit out¬ 
side of the “postponement concept. Firms in the food 
sector dynamically move from a MTS to MTO status 
to manage short term changes in demand operating a 
hybrid MTS/MTO approach not based on “postpone¬ 
ment”. 

MTS/MTO processes enjoy popularity in industry 
and find themselves implemented in software sys¬ 
tems. Enterprise resource planning systems such as 
SAP, Microsoft Dynamics and Oracle’s JD Edwards 
accommodate strategies for MTS and MTO prod¬ 
ucts. Eor instance SAP has a detailed default process 
flow for handling MTO products. The user is forced 
to distinguish a product to be either of MTO or MTS 
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type. That means products of hybrid character lead to 
two product records, which evolve into two demand 
and production time series. In this paper we develop 
an alternative approach that reverses the process by 
having one demand time series for a hybrid product 
and splitting it up into MTO and MTS items for pro¬ 
duction planning and control. Some researchers have 
begun tackling the challenge of modeling MTS/MTO 
systems and to analytically decide whether a prod¬ 
uct should be classified as MTS or MTO through the 
lens of demand analysis. The methodology section 
builds on this earlier approach through the applica¬ 
tion of outlier detection techniques leading to the in¬ 
troduction of the new method developed to dynami¬ 
cally identify MTO and MTS product categories. 


3.2. Outlier detection techniques 

Demand analysis has been utilised to cate- 
gorise products into MTS an d MTO families. 
D’Alessandro and Baveja ( 2000l) suggested that an 
analysis through the prism of demand variability and 
average weekly demand products can be categorised 
as MTO or MTS. A fair amount of literature in this 
field utilises the coefficient of variation (CoV): 


cr 


( 1 ) 


which is the ratio of standard deviation cr and mean 
/i. Products with low demand variability and “high 
volume” were classed as MTS, remaining products 
were viewed as MTO. A typical characterisation of 
a MTO product is its low average demand and high 
variation.Some researchers have stated that a coef¬ 
ficient of variation on or below 0.5 may be used 
to classify a product as MTS. The authors have not 
found any definitive method in the literature, which 
computes the likeliness of a product falling into the 
MTS or MTO category. Determining where the cut¬ 
off point is for demand variability and the portion¬ 
ing of products into MTO or MTS based on the 
CoV factor can be subje ctive and difficult to deter¬ 
mine (ISoman et al.Ll2007h . Recent research has chal¬ 
lenged the validity of applying of CoV as a lens 
to allocate MTS or MTO status due to the frequent 
changes in CoV for products as market and their de- 
man ds change quickly altering a prod uct categoriza¬ 
tion ( Godsell and Kharlmanovl 2012 ). The low av¬ 


erage demand could be deceiving because large, un¬ 
forecasted, orders can shift the average to any level. 
The ratio of the sample standard deviation to the 
sample mean overcomes this problem. However, 
trends can affect this ratio. Removing the linear trend 
leads to improvements. 

In practice operations need to derive a production 
strategy that accommodates MTS and MTO from a 
single demand profile. We will assume that outliers 
represent MTO items. Henc e, we have to consid¬ 


ered outliers in time series. iLewisI (119751) reviews 


demand analysis in his extensive work where he dis¬ 
cusses demand impulses in the context of adaptive 
forecasting He emphasises that “unusual demands 
still have to be deal t with’" indicating the relev ance 
of outlier detection. ISilver and PetersonI (119851) also 
acknowledge that “temporary change in demand pat¬ 
tern"' such as promotions should be filtered out be¬ 
fore forecasting.Removing identified outliers will re¬ 
duce the variation of the MTS classified product thus 
re ducing subsequ e nt inv entory. 


Atkinson et al. ( 1997b suggested ways of check¬ 


ing fitted models to possible shocks and introduced 
a new type of intervention analysis, which uses an 
extra parameter for the outlier. This is equivalent to 
deleting the outlier and observing the residual sum 
of squares in a regression model. A considerable 
a mount of lite rature about outlier detections refers 
to IFoxI (1972), which introduces likelihood ratios. 
Fox looks at stationary time series and proposes two 
types of errors. The first type deals with recording 
error or “isolated independent cross execution” that 
are independent of other observations. The second 
type of “outliers” affects consecutive observations - 
known as additives or innovations. Our work focuses 
on the first type, which is commonly addressed as 
shock in the time series literature. Time series analy¬ 
sis usually relies on predictive causality, the opposite 
o f what we aim to do. 

Hawkins ( 198fll) discusses the identification of out¬ 
liers in great detail and defines them as: “An outlier 
is an observation which deviates so much from the 
other observations as to arouse suspicions that it was 
generated by a different mechanism”. It is conve¬ 
nient to introduce a basic model, which is univariate, 
assumes one “normal” generatin g mechanism and 


that outliers are rare observations. Barnett and Lewis 
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(119941) similarly define an outlier in a set of data to 
be “an observation (or subset of observations) which 
appears to be inconsistent with the remainder of that 
set ofdata’\ On the other hand contaminants are ob¬ 
servations which are not “genuine members” of the 
main distribution. Extremes could be outliers or con¬ 
taminants. Assume that the actual distribution F is 
known for a sorted sample xi, X 2 ,..., Xi and 
are the sample extremes, which might be outliers. 
However, outliers are extreme values. Contaminants 
interfere with the original distribution F and may 
stem from another distribution G. This could result 
in extreme values wrongly identified as outliers. An 
intuitive model could assume that a demand time se¬ 
ries is made up of a MTS distribution F and a MTO 
distribution G. 

A variety of methods have been devised to de¬ 
tect outliers. Typically these techniques fall into 
two categor ies: model based an d proximity based 
approaches ( Kriegel et al.l 2010l) . Model based ap¬ 
proaches can use statistical tests, depth considera¬ 
tions and deviation analysis. Proximity based tech¬ 
niques usually use distan ce measures or density as¬ 
pects. Chen et al. ( 2010h compare outlier detection 
techniques in regards of of performance. They anal¬ 
ysed statistics-based, distance-based and density- 
based approaches. Model based approach e s hav e 
been extensively investigated by HawkinsI (Il980b . 
Most of these models assume normal distributed 
data; however, all samples in the case study re¬ 
ject this assumption. We will focus our review on 
proximity based techniques. Distance based ap¬ 
proaches have the basic assumption that data has a 
dense neighbourhood. This allows identification of 
outliers due to their distance frorn the neighbour¬ 
hood. Knorr and Ngl ( 1998 , 1999 ) have discussed 
such algorithms based on iodizes, nested-loops and 
grids. The k-nearest neig hbour method was adapte d 
for outlier detection by Ramaswamy et al. ( 2()0()h . 
Angiulli and Pizzutil (120021) and many more. Prox¬ 
imity techniques based on density usually assume 
that the densit y around an outlier is diff erent to 
its neighbours. Breunig et al. ( 19991 2000h use the 
concept of the Local Outlier Factor that overcomes 
the issue of clusters having different den sities. Re¬ 


lated t o this concept a re the works from iTang et al. 


(2002) and Jin et al.l ( 200 lb . A further improve¬ 


ment was achieved by I.Jin et al.l (120061) by using 


symm etric neighbourhood relationships. ICao et al. 


( 2010b developed a density-based algorithm to iden 
tify outliers. They introduced a density-similarity- 
neighbourhood factor which suggests the likeliness 
of the data to be an outlier. Another interesting ap¬ 
proach which uses a local outl ier correlation inte¬ 


gral o n an e neighbourhood is iPapadimitriou et al. 
( 2003 ) work. Kriegel et al. ( 2010b mentioned that 
cluster algorithms can be used to identify a noise 
set/outliers. Multiple outliers similar to each other 
would form a cluster r ather than being identified as 


noise. 


Alan and Catall (2011) use thresholds to form 


multiple clusters using proximity based outlier de¬ 
tection mechanisms. These techniques are insightful 
but were not designed for demand time series in par¬ 
ticular. 


3.3. MTO cluster detection technique 

In this work we propose a new technique which 
is based on proximity and density. The novelty of 
our approach is a distance measure operating on the 
probability mass distribution that identifies MTS and 
MTO categories via clustering. Before we embark 
on the details of this technique we give a first demand 
taxonomy of hybrid MTS-MTO products: 

1. MTO outliers 

(a) Demand is high and infrequent 

(b) Demand is low and infrequent 

2. MTO contaminants 

(a) MTO and MTS demand distributions dif¬ 
fer 

(b) MTO and MTS demand are similar 

Type 1 (a) assumes that MTO items are characterised 
by high demand and low frequency. This is the case 
we will develop our techniques for. Type 1 (b) is 
not important to us because we assume that MTS 
production will not be affected significantly as low 
demand can be served from inventory. Type 2 (a) 
allows a separation of the distributions when deter¬ 
mination of the parameters is a possibility. Type 2 
(b) constitutes a challenge and cannot be analysed 
without further information and will be discussed in 
future research. 
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4. Methodology 


In this section we will introduce new methods, 
which help to classify what parts of the demand 
should be used as MTS or MTO. Some of the the¬ 
oretical underpinnings were introduced in the previ¬ 
ous section, which motivated a discussion of outlier 
detection techniques. Furthermore we derived a tax¬ 
onomy of hybrid MTS-MTO systems. 

Here, we will propose our new MTS and MTO 
identification method. This requires us to introduce 
a novel series distance measure, which operates on 
histograms and time series. We embed this distance 
measure into a cluster technique and derive three cat¬ 
egories. A new immediate neighbour detection tech¬ 
nique reduces these to the two final MTS and MTO 
categories. 

The methodology introduced in this research fo¬ 
cuses on type 1 demand profile. Our goal is to sep¬ 
arate the convoluted MTS-MTO demand. The first 
step is the definition of an appropriate distance met¬ 
ric. A new city block metric that operates on a his¬ 
togram with frequencies / = (/i, ...,/„) is proposed. 

dst = \fs - ft\ + Q:\Cs - Ct\, (2) 


The distance d^t between s and t uses equidistant 
classes c = (ci,..., c„) with d = c,+i - c,- > 0,i < n, 
where Cj represents the centre of class i. Note that 
a scaling factor a was introduced. The factor is sup¬ 
ported by the following considerations. Assume each 
class contains only one observation then the unit dis¬ 
tance between classes seems appropriate. If each 
class contains ten observations a unit distance would 
be inappropriate, because the class distance is not in 
proportion to the frequency scale. This motivates the 
class scaling factor a: 


ZLi ft 

(n - l)d 


( 3 ) 


Example 4.1 (scaled city block distance). Let fig¬ 
ure 121(a) be the histogram derived from a demand 
time series. Class two shows demand for two items 
occurs 80 times/days. A total of 195 days were ob¬ 
served and placed into seven demand classes. The 
scaling factor is a = 195/(7 - 1) = 32.5. Hence 


the distance between s = (2,80) and t = (5,20) is 
dst = |80 - 20| + 32.512 - 5| = 60 + 97.5 = 157.5. 
Note that the demand classes can be transformed to 
represent other quantities (e.g. c\ could represent 
1200 items, ci = 1300 with d = 100)). 


The problem with the above distance measure is that 
for instance the distance between two classes sepa¬ 
rated by another one could be shorter than the one 
next to each other. As a result clusters could consist 
out of disconnected classes. We have overcome this 
by introducing a novel “series” distance measure. 
The distance between classes next to each other is de¬ 
termined as before, i.e. dk,k+i = \fk-fk+i\+Q:\Ck-Ck+il 
Note that a\ck - Ck+i \ can be simplified to the average 
frequency /. The distance between two classes s and 
t with 5 < t is determined by adding up previous dis¬ 
tances: 

t-i 

dst = '^^dk,k+\- (4) 

k=s 


This can be expressed as a recursion dk,k +2 = dk,k+i + 
dk+i,k+ 2 - All individual series distances form a corre¬ 
sponding matrix D = {dsf. Example |2ls histogram 
has similar characteristics to those observed in the 
case study, i.e. classes with low demand occur more 
frequently than those with high demand (a type 1 
MTS-MTO hybrid). These distance measures are 
used to build an agglomerative hierarchical cluster 
tree. We use the Unweighted Pair Group Method 
with Arithmetic M ean (UPGMA) to bu i ld a d endo- 
gram (binary tree) ( Sokal and Michen^ 1958 ). 


Example 4.2 (Dendogram). We continue the previ¬ 
ous example and build an agglomerative hierarchical 
cluster tree (see figure^(b)). The average distance 
between class 1 and 2 is 62.5. The average distance 
to class three from one and two is = 73.8. This 
cluster represents MTS items. Class 5 is in between 
the MTS and MTO cluster. 


In general the distance of two clusters S and O are 
determined by the average pair distances: 


1 

nm 


ZZ4,. 

seS oeO 


( 5 ) 


where n and m are the number of elements in S 
and O respectively. A possibility to improve the 
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Figure 2: Example (a) demand histogram showing scaled city block metric, (b) dendogram. 


distance measures further is by eonsidering holding 
eosts, setup eosts and shortage eosts. 

The proeedure developed above is summarised in 
algorithm [TJ This proeedure offers a systematie way 


Algorithm 1 MTS/MTO eluster deteetion 
Require: equidistant time series of demands y, no 
major innovation in time series 
Ensure: MTO threshold m, MTS threshold M 

1; remove linear trend from time series y and obtain 

2: transform y into histogram with elasses C; and 
eorresponding frequeneies / 

3; ereate series distanee matrix D = (dst) using 
equation Q 

4: determine hierarehieal eluster tree using un¬ 
weighted average distanees (see example 1431) . 

5: extraet three elusters from the tree, identifying 
MTO and MTS thresholds 


of identifying initial MTS and MTO thresholds. It 
may be the ease that a produet has no MTO produets. 
However, the above proeedure will always identify a 
potential MTO elass. In sueh eases we suggest to 
fall baek to elassie serviee level strategies. For in- 
stanee ensuring that predieted demand is fulfilled in 
95% of the time periods. We require that there are no 
“major” innovations in the time series. However, that 
issue and a solution are addressed in the following 
seetion. The first step is the remove of linear trend. 
This is a first approximation whieh ean be followed 


with higher order trend removals in aeeordanee to a 
Taylor series development. Next a histogram is de¬ 
rived. We would like to sugges t as future work to 
use Kernel Density Funetions. Botev et al. ( 2010h 
have proposed a method whieh has the advantage of 
being free of an underlying normal distribution as¬ 
sumption. Step 3 and 4 were already explained in 
this seetion. The last step identifies the MTO and 
MTS thresholds by using the agglomerative hierar¬ 
ehieal elustering that identifies three elusters. The 
introdueed distanee ensures that the elusters are ad- 
jaeent rather than overlapping. 


4.1. Time series 

In the previous seetion we have introdueed a teeh- 
nique to identify MTO items; assuming a hybrid 
MTS-MTO time series. Removing MTO items from 
the hybrid time series will improve subsequent time 
series analysis for the remaining MTS items. This in 
turn eauses savings in the holding eosts. The produe- 
tion sehedule ean be determined via solving the Eeo- 
nomie Lot Seheduling Problem. Alternatively gen¬ 
eral time series analysis and foreeasting teehniques 
ean be used to determine the antieipated demand and 
derive a produetion sehedule. 

Let us eonsider three theoretieal eases of time se¬ 
ries shapes: (1) eonstant demand; (2) linear inereas- 
ing demand; and (3) demand undergoes an innova¬ 
tion. The first ease leads to histograms with a sin¬ 
gle elass of high frequeney, whieh is identified as 
MTO or MTS. This is in aeeordanee to our empir- 
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ical expectations, i.e. either orders are known in 
advance or excellent predictability. That means no 
further information is required to make a decision. 
The second case leads automatically to case one af¬ 
ter trend removal. Not removing the trend would 
equipartition the histogram, suggesting wrongly that 
higher demand consists out of MTO items. The third 
case is a change in trend (innovation), which linear 
trend removal would not cope with. However, an 
additional preprocessing step to the previously intro¬ 
duced procedure [His able to transform case three into 
case one or two. Finding the innovation is achieved 
by using a weighted average intra-cluster distance 
measure (weighted pair group method with averag¬ 
ing, WPGMA) instead of the average measure. First 
the series distance (equation (4]) is determined for the 
time series y. This is followed by the WPGMA clus¬ 
tering and identifying two clusters. 

The boundary between the two clusters identifies 
the innovation. This is the cut-off point, that means 
we can use the remaining observations for cluster¬ 
ing and identifying MTO items using algorithm [T] 
Hence, we will improve the MTS demand predic¬ 
tions for the shortened time-series. Certainly the 
question remains open, what is a “major break”? Our 
method divides the time series into two parts, which 
allows several approaches to find an answer. Two 
approaches could be used to answer the question: 
( 1 ) comparing the angle between the trend lines; ( 2 ) 
expectation and variance considerations of the two 
parts. 

4.2. Summary 

An decision support system was developed that 
classifies MTO and MTS items. To be more spe¬ 
cific the proposed method transforms a time into a 
histogram. The histogram is further processed us¬ 
ing a cluster detection technique that identified MTO 
and MTS categories. The success of this clustering 
is ensured by the proposed distance measure. 

In practice there are several issues, which still 
need to be dealt with. The next section will look at 
those and add a few additionally required techniques, 
which are better understood through the application 
of actual data. 


5. Case study results 

The background of the case study company was 
given in the introduction and elaborated on the fact 
that there are 250 products to be manufactured. Ini¬ 
tially we will focus on one particular product, and 
later discuss a set of ten products which were pro¬ 
duced on the same production line. This will make 
use of the methodology introduced in the previous 
section. This section will give us a good understand¬ 
ing of the data, and introduce one more required 
technique motivated by the case study’s practical re¬ 
quirements. This technique is an approach for prac¬ 
titioners and suggests an approach to dissolve the 
MTS or MTO “grey zone”. 

5.1. Product categories 

Figure [3] (a) shows the demand profile for prod¬ 
uct GOISL over a period of two years (n = 104 
weeks). It displays the frequencies of MTS/MTO 
items in 20 classes (bins). This histogram was de¬ 
rived from the time series displayed in figure |4] (a). 
The time series consists of equidistant demand ob¬ 
servations y = yi,y 2 , • • • i-e. y,- is the number of 
items moved from the central warehouse (see section 
O to the customer in week i. A high variability can 
be observed. A moving average of four weeks (black 
line) is shown and approximates a monthly demand 
profile. The linear trend (red line) shows that the de¬ 
mand is non-stationary. 

The coefficient of variation ^ = 1.1 suggests a 
product being predominantly MTS. Here x is the 
sample mean and s is the sample standard devia¬ 
tion. MTO items are visible by computing the z- 
score z = ^. This identifies two weeks as outliers 
which have a z-score above 3.1. They are highlighted 
in figure |4] (see dashed red lines). To be more pre¬ 
cise week 73 and 82 have a z-scores of Z 73 = 6.5 
and Z 82 = 3.8 respectively. That means for a normal 
distribution this would not be expected in 99.9% of 
the cases. Despite the correct identification it should 
be noted that the KolmogorovSmimov test does not 
confirm a normal distribution. 

Considering figure [3] the demand can be divided 
into three categories: MTO, MTS and a third one. 
Algorithm [T] identifies these categories using the hi¬ 
erarchical agglomerative clustering (see section (33]). 
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(a) (b) 


Figure 3; Product GOISL (a) Probability density and mass of demand identifying MTO and MTS, (b) Dendogram. 


The clusters were created using the series city block 
distance (equation ID). Figure [3] (b) visualises and 
quantifies the intra cluster distances. The MTO cat¬ 
egory occurs in two weeks with quantities 25,650 
(12.8%). MTS comprises 88.5% weeks and a to¬ 
tal quantity of 125,041 (62.3%). The third category 
(MTO or MTS) happens in 9.6% weeks with a quan¬ 
tity of 49,935 (24.9%). 

Our objective is to get clarity about the third cate¬ 
gory. In order to achieve this we have to analyse the 
individual weeks. We will use the symbol and 
Wo to denote the MTS and MTO weeks respectively. 
The union of all weeks is abbreviated with W. That 
means the weeks within the MTS or MTO category 
are W \ {W„ U Wo). The weeks can be derived using 
the MTO and MTS demand threshold determined by 
algorithm [H Ultimately the intention is to end up 
with the MTS and MTO categories only. We can as¬ 
sume that the MTS category is pure with an average 
demand of 1,359 items per week. The two weeks 
identified as MTO occur in week 73 and 82 with de¬ 
mands of 15,700 and 9,950 respectively. A closer 
investigation of the MTO neighbours helps in decid¬ 
ing about which of the weeks W)j should be added 
to the current MTO weeks Wo- Week 73 has only 
one neighbouring week with unusual high demand 
5,096 (classified within the MTS or MTO category). 
The other neighbouring weeks [69,76] \ {73,74} have 
an average demand of 2,899. A similar scenario is 
happening with week 82. Week 81 has a demand 
of 6,663 items. The surrounding 7 weeks window 
[78,85] \ [81,82] has an average demand of 1,542. 
The immediate neighbours of MTO from the MTS 
or MTO category have a demand volume of 5.9%. 


We conclude that close to MTO identified weeks are 
likely to be part of MTO category. The pure MTS 
sample mean underestimates the MTO neighbouring 
demand and cannot be used for estimating the de¬ 
mand during the MTO weeks. These heuristic con¬ 
siderations are used to decide about the final MTS 
and MTO classes. 


5.2. Final classification technique 

Based on the above considerations we will intro¬ 
duce a decision support mechanism to get even more 
clarity about MTS and MTO demand. We begin by 
using the third category (containing MTS or MTO 
weeks) as candidates for the MTO category. A can¬ 
didate ai is the tuple (w,-,y/), where w; is the week 
and y,- the corresponding demand. If a candidate is 
in the immediate neighbourhood (i.e.w,_i or W/+i) of 
a MTO identified week then the candidate becomes 
part of the MTO category. This is followed up recur¬ 
sively, i.e. a newly assigned MTO (ex-candidate) can 
cause candidates to be added to the MTO category. 
The set of MTO-candidate weeks identified this way 
will be abbreviated with Wf For instance for prod¬ 
uct GOISL we only have two candidates, which are 
added to the MTO category. Now we are in the posi¬ 
tion to create our final demand profiles for MTO and 
MTS items. The final MTO weeks are Wo := WoUW^ 
and the MTS weeks are := W \ Wo. We assume 
that weeks with MTO demand also have MTS de¬ 
mand, which is equal to the average MTS demand 
Ds. Replacing the demand of the MTO weeks Wo 
with Ds gives us the time series displayed in figure SI 
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(a) 


(b) 


Figure 4; Product GOISL (a) hybrid time series, (b) MTS time series due to removal of MTO items. 


which we abbreviate with y^. 


yl 


\d„ ifieWo 


( 6 ) 


The other time series is obtained by taking the de¬ 
mand during the MTO weeks and reduce it by 
Ds. It is zero except for the MTO weeks and its re¬ 
cursively added eandidates. 


0 




if ieW\Wa 


\yi-D„ if ie Wo 


(7) 


Of eourse -l- y° is equal to y,. 

Considering product GOISL the final MTS de¬ 
mand is 84.6%, which is an increase of 22.3% to 
the pure MTS cluster demand. The MTO compo¬ 
nent has an apparent inerease by 2.6% resulting in 
a demand of 15.4%. We use the word “apparent” 
because the MTO cluster demand eontains demand, 
whieh we have estimated to be per week. 

In summary this additional assignment teehnique 
has created a time series for MTS units and one for 
MTO units. 

5.3. Practitioner approach and ten products 

We will give a guideline for praetitioners to deeide 
about MTO and MTS products/items in the absenee 
of programming tool. These guidelines were applied 
to ten produets. The timeseries for the ten produets 
produeed on the same manufacturing line are dis¬ 
played in figure [51 Results in regards to savings and 
CoV are diseussed in this eontext. Here savings are 


due to reduced holding cost, because the production 
knows the fraetion of MTO units. The argument is 
that MTO items can be produced just in time. That 
means an improved planning process is possible. 

The approaeh for practitioners is shown in algo¬ 
rithmic! As a first step the CoV should be used. Table 

Algorithm 2 Approach for practitioners 
Require: Time series of demands D 
Ensure: MTO threshold m, MTS threshold M 
1; Analyse time series 

2: if no break in structure (innovation) “visible” 

then 

3: Identify MTO and MTS thresholds by 

4; Create histogram with frequeneies / and 
quantity eentres c,. 

5: Determine deltas A between eonseeutive fre¬ 

quencies - / 2 :„. 

6: MTS threshold M is the quantity after the 

maximal frequency change max A 
7; MTO threshold m is the quantity after the 
next maximal frequency change 
8; else if break in structure identified then 
9; additional analysis required (set traffie light 
to orange) 

10; else if eoefficient of eovariation high then 

11: MTO likely (alternatives: new product, 

product discontinuation) set traffie light to red 

12 : end if 

iHgives the CoV for the ten produets. Two of them re¬ 
veal a CoV significantly above the average (product 
4 and 10). The CoV was determined following the 
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Figure 5; Time series showing the weekly demand. 


removal of the linear trend as discussed in section 
4.1. Typically the removal of trend reduces the CoV, 
that means without removal a identification of MTO 
products is less reliable. Closer investigation of the 
two products reveals one as a new product being in¬ 
troduced and the other one as a pure MTO product 
with low demand. We will identify them with a red 
traffic light, in table 1, to indicate that no further in¬ 
vestigation into these two products will take place 
(see algorithm[2l third if clause). 

As a second step the practitioner can visually iden¬ 
tify those time series showing an obvious innovation 
such as a break in structure or change in trend. These 
are marked with an orange traffic light in table [B In 
our case study product six, eight and nine were iden¬ 
tified. Algorithm |2] only suggests that further anal¬ 
ysis is required. In subsection 14.11 we introduced a 
time series clustering technique that allows to find 
the innovation for such time series. Using the new 
method allows the remainder of the time series to be 
analysed. Product six and eight remaining demand 
profiles were identified as hybrid MTS-MTO. Their 
estimated savings due to separate MTS and MTO 


production are 8.9% and 15.0%. Product nine is clas¬ 
sified a MTO and no savings are possible. 

The third step deals with the hybrid MTS-MTO 
character of the demand profiles. It should be noted 
that the proposed procedure omits or changes steps 
addressed in the original algorithm [TJ The major 
changes are the trend removal and the clustering 
technique. These adaptations lead to differences in 
the savings of about 5%. The results obtained with 
algorithm[T]are shown in table [B and [2l As we can see 
the savings are 18.1% on average. This is due to the 
reduced holding costs by assuming neglect able hold¬ 
ing costs for MTO items. That means MTO items are 
assumed to be produced just-in-time. Figure visu¬ 
alises the MTS/MTO categories. The last sub-plot 
displays the coefficients of variation and their group¬ 
ings. The first table shows that the estimated savings 
are in the range of 13.4% and 27.0%. Furthermore 
the CoV reduced in all cases (compare CoV and 
CoV-S), which will result in additional savings due to 
the reduced variation and increased production sta¬ 
bility. Here, CoV-S is the abbreviation for the coeffi¬ 
cient of variation for the MTS products, which varies 
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Figure 6; histograms demonstrating the MTS, MTS or MTO and MTO identihcation (trend removal when green traffic light). 


Table 1: Savings for products due to MTO item identihcation. 


Product 

TL 

Savings 

#weeks 

CoV 

CoV-S 

CoV-O 

BO 

Z 

GOISL 

• 

15.4% 

104 

1.10 

0.79 

5.70 

0.04% 

2 

G02SN 

• 

13.4% 

104 

0.63 

0.50 

3.38 

0.00% 

2 

G03MH 

• 

22.6% 

104 

1.20 

1.07 

3.83 

1.65% 

2 

R04TN 

• 


104 

5.13 


5.13 

0.00% 

4 

G05MM 

• 

27.0% 

104 

0.98 

0.77 

2.96 

0.09% 

2 

O06MN 

• 

8.9% 

44 

0.52 

0.45 

3.66 

0.00% 

0 

G07ML 

• 

24.2% 

104 

0.93 

0.77 

2.96 

0.02% 

2 

O08SN 

• 

15.0% 

48 

0.40 

0.32 

1.98 

0.00% 

0 

O09LN 

• 


43 

1.73 

4.86 

1.73 

0.00% 

0 

RIOLN 

• 


104 

3.86 


3.86 

0.00% 

4 


Average 

18.1% 

86.3 

1.65 

1.19 

3.52 

0.18% 

1.8 


between .77 and 1.07. CoV-0 is the CoV for MTO 
produets and varies between 2.96 and 5.70. A teeh- 
nieality should be noted; there are weeks with neg¬ 
ative demand (baek-orders, BO). In particular prod¬ 
uct three had a negative demand of 1.65%. For the 
clustering technique to work effectively back-orders 
above 0.3% were treated separately. We have dis¬ 
played the number of weeks as outlier via the z-score 
method in the last column of table [B These numbers 
are equal or below the MTO weeks which were iden¬ 
tified via algorithm [TJ The z-score assumes a nor¬ 
mal distribution; however, all ten products failed the 
Kolmogorov-Smimov and the Jarque-Bera test for 


normality. Many outlier detection methods depend 
on the normal distribution (see section 1X21) . The em¬ 
pirical data cannot support their relevance for this 
case study. This is in accordance with time series 
literature in general, which attempts to find as a last 
remainder normally distributed noise. Table [2] shows 
the actual number of items per category. 

Clearly indicating that product one and two have 
significantly more demand than the other products. 
The table shows first the absolute total demand over 
two years (except for the products with orange traffic 
lights). These numbers were obtained using algo¬ 
rithm dl They are also represented as volume per¬ 
centage. Using the recursive immediate neighbour 
approach the MTS or MTO category was dissolved 
resulting in the “final” two MTS and MTO cate¬ 
gories. An average of 81.9% MTS items were found 
in the hybrid products. The last column of this table 
shows the average absolute number of MTS items in 
the final MTS category. 

5.4. Comments 

The food processing sector continues to face in¬ 
creasing product volatility and service heterogene- 
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Table 2: MTS and MTO volume per category. 


Product 

TL 

MTS 

MT-S/0 

MTO 

MTS% 

MT-S/0% 

MTO% 

fMTS 

fMTO 

fMTS% 

fMTO% 

mMTS 

GOISL 

• 

125,041 

49,935 

25,650 

62.3% 

24.9% 

12.8% 

169,746 

30,880 

84.6% 

15.4% 

1,632 

G02SN 

• 

183,955 

145,144 

52,282 

48.2% 

38.1% 

13.7% 

330,254 

51,127 

86.6% 

13.4% 

3,176 

G03MH 

• 

12,088 

17,737 

8,416 

31.6% 

46.4% 

22.0% 

29,582 

8,660 

77.4% 

22.6% 

284 

R04TN 

• 



320 



100.0% 


320 


100.0% 


G05MM 

• 

19,246 

22,581 

17,401 

32.5% 

38.1% 

29.4% 

43,224 

16,004 

73.0% 

27.0% 

416 

O06MN 

• 

15,714 

18,752 

4,202 

40.6% 

48.5% 

10.9% 

35,217 

3,452 

91.1% 

8.9% 

800 

G07ML 

• 

13,161 

14,695 

10,573 

34.2% 

38.2% 

27.5% 

29,138 

9,291 

75.8% 

24.2% 

280 

O08SN 

• 

35,463 

96,268 

23,550 

22.8% 

62.0% 

15.2% 

132,044 

23,237 

85.0% 

15.0% 

2,751 

O09LN 

• 



14,150 



100.0% 


14,150 


100.0% 


RIOLN 

• 



6,576 



100.0% 


6,576 


100.0% 



Average 

57,810 

52,159 

16,312 

38.9% 

42.3% 

43.2% 

109,886 

16,370 

81.9% 

42.7% 

1,334 


ity altering the historieal demand profiles of SKU’s. 
Seheduling and eontrolling produets on the basis of 
a pure MTS approaeh has found in the ease study 
firm to be futile and eostly. Historieally eategorised 
MTS produets regularly exhibited eharaeteristies that 
align with a MTO approaeh rendering the produc¬ 
tion schedule redundant. Rush and exceptionally 
large orders were becoming the norm for the case 
study firm and were found to be difficult to manage 
through standard planning and control procedures. 
The methodology outlined in this section represents 
the modeling, analysis and subsequent implementa¬ 
tion of a new way of planning and controlling prod¬ 
ucts that demonstrate duality in their product cat¬ 
egorisation. The algorithm offers practitioners the 
possibility to evaluate products and demand in terms 
of the production approach to be utilised to satisfy 
demand. The methodology above provides a support 
tool for decision-making and scheduling of produc¬ 
tion at the SKU level that separates orders in terms of 
MTO and MTS .The result for food processing firm 
showed that a separate production of the MTO items 
led to a reduction of holding costs by 18.1% on aver¬ 
age. 

6. Discussion and conclusion 

This empirical research has demonstrated how to 
deal with the combined MTS-MTO situation in the 
context of finite capacity. To date the authors are 
not aware of the existence of a published quantita¬ 
tive method which identifies MTO and MTS items 
automatically within a hybrid demand time series. 
The paper has reported on the challenges and diffi¬ 
culties of splitting MTO and MTS demand profiles 


to enable firms to manage fluctuations in customers 
demand and profile. The developed methodology 
provides practitioners with the tools to revise their 
order history and correlate it with the demand time 
series. Specifically the method supports the identi¬ 
fication and scheduling of change in product status 
between MTS, MTO and MTO/MTS that could im¬ 
prove customer service and reduce costs. The new 
methodology has the potential to improve the per¬ 
formance of production planning and control within 
the food processing sector. The algorithm developed 
for practitioners is the first step in aiding production 
planners in achieving a stable schedule when faced 
with products exhibiting shifting categorisation sta¬ 
tus. 

The developed methodology is a valuable contri¬ 
bution to the MTO-MTS literature as it addresses, 
empirically as well as mathematically, a problem that 
planners face on a regular basis. The interaction ef¬ 
fects of MTS-MTO demand characteristics with fi¬ 
nite capacity is an area that had been identified by 
several researchers as an important area for research 
in planning and control. The data analysis from the 
study of ten products paints a complex picture of the 
dualistic nature of products in the food processing 
sector and suggests that further research is required 
to develop a stronger link between the developed 
methodology and practice. Identifying the MTO or 
MTS status of demand and converting the informa¬ 
tion to develop a robust manufacturing schedule, in 
the context of a dynamic MTO/MTS environment, is 
an area for future investigation. 

Through the development and implementation of 
the new methodology two further research directions 
have emerged. MTO and MTS demands that are sim- 
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ilar in distribution are a challenge that have still to 
be researched and resolved in terms of the model. 
Identifying, segregating and managing the overlap in 
distributions requires further data collection over a 
wider range of products and longer time-series. Inte¬ 
grating the single demand series for a hybrid product 
into the enterprise resource planning (ERP) system is 
also highlighted as an area for future research. Given 
the wide application of ERP systems in the food pro¬ 
cessing sector linking the algorithm to the main tool 
of production planners is a vital next step in develop¬ 
ing a pathway for successful application and imple¬ 
mentation of the new methodology. 
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